Fast and slow edges in bilayer graphene nanoribbons: 
Tuning the transition from band- to Mott-insulator 
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We show that gated bilayer graphene zigzag ribbons possess a fast and a slow edge, characterized 
by edge state velocities that differ due to non-negligible next-nearest-neighbor hopping elements. By 
applying bosonization and renormalization group methods, we find that the slow edge can acquire 
a sizable interaction-induced gap, which is tunable via an external gate voltage V g . In contrast to 
the gate-induced gap in the bulk, the interaction-induced gap depends non-monotonously on the 
on-site potential V. 
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I. INTRODUCTION 



One of the most attractive properties of bilayer 
graphene (which is made of two coupled atomic layers 
of carbon) is the ability to induce a tunable spectral gap 
Ay by applying a perpendicular electric field»i~— Rapid 
advances in patterning graphene on the nanoscale now 
make it feasible to fabricate graphene ribbons with well 
defined edge termination^^ while experiments on nar- 
row ribbons show that they can display a gap Aw due 
to transverse size quantization^ The ensuing facility to 
confine electrons in a controllable way via gate poten- 
tials and patterning makes bilayer nanoribbons a promis- 
ing candidate for nanoelectronic applications. Systems 
suitable of being used as a basis for a transistor should 
exhibit a high resistance in the off state. In the case 
of nanoribbons with zigzag termination, however, this 
requirement poses a problem, since the edges support 
current-carrying states with energies inside the bulk gap, 
and a much reduced hybridization gap A/, <C Ay, Aw 
due to exponentially small tunneling between the edges 
2r— On the other hand, several novel concepts (notably, 
valleytronicsi^) exploit the existence and specific proper- 
ties of such edge states. 

Previous works on edge states in gated bilayer ribbons 
ignore two elements: (i) next-nearest-neighbor hopping, 
which is known to strongly affect the bulk properties^— 
and is of crucial importance for edge states in mono- 
layer ribbons^ and (ii) the effects of interactions, which 
are known to strongly influence the properties of quasi- 
one dimensional systems, including carbon-based sys- 
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In this paper, we point out that next-nearest- 
neighbor hopping breaks the symmetry between the 
edges, resulting into a slow edge and a fast edge charac- 
terized by different values of the propagation velocity. We 
then explore the consequences for the question whether 
many-body effects can help to open a gap in these chan- 
nels. Treating interactions in a Hubbard model on the 
basis of bosonization and renormalization group methods 
at T = 0, we find that the edge channels display a Mott 
transition at half filling. In presence of the higher-order 
hopping terms, the associated charge gap A p on the slow 



edge can take on sizeable values, while on the fast edge 
the gap is negligibly small. Remarkably, the gap depends 
sensitively on the applied gate voltage V; the effects of 
interactions can therefore be controlled externally. Fur- 
thermore, the dependence of A p on V is non-monotonic, 
in striking contrast to the behavior of the field-induced 
gap Ay in the bulk of the system. 

This paper is organized as follows. In Sec. II we 
present the hamiltonian for the graphene bilayer rib- 
bon including the next-nearest-neighbor hopping and the 
hopping term between non-dimerised atoms in both lay- 
ers. We then describe the edge states for a zigzag rib- 
bon of bilayer graphene and their different localization 
properties along the transverse direction of the ribbon by 
computing the Inverse Participation Ratio. In Sec. Ill 
we describe the effects of considering electron-electron in- 
teractions in the spectrum for edge states. We find that 
because the different properties of the two edges the gap 
opened by the interaction is different for both edges. We 
also briefly discuss a possible experimental verification of 
this different behavior. We summarize our results in Sec. 
IV. 



II. TIGHT BINDING MODEL 

A graphene bilayer nanoribbons with zigzag edge ter- 
mination is shown in Fig. [T] (left panel). Focussing 
on the low-energy bands which participate in electronic 
transport, such ribbons can be modelled by a tight- 
binding hamiltonian H = E<T=jf Eij7ijV c J>i where 
each carbon atom hosts one spin-degenerate electronic or- 
bital with annihilation operator Ci a . Recent experiments 
specifically addressed the values of the parameters jij in 
bilayer graphene J^r— Nearest-neighbor hopping in the 
same graphene layer is described by a hopping element 
—70, where 70 — 3eV. The layers are Bernal stacked, 
with hopping element 71 ~ 0.1 270 between dimerised 
atoms that lie on top of each other. The symmetry be- 
tween both layers can be broken by top- or back-gates 
with a gate voltage V g , which induce a perpendicular 
electric field and give rise to an on-site potential ±V on 
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FIG. 1: (Color online) Upper panel: Sketch of a bilayer 
graphene nanoribbon of lattice constant a and width W (in- 
set: definition of hopping terms). Middle panel: Dispersion 
relation of edge states for a nanoribbon with W = 32 x y/3a 
and layer-symmetry breaking potential V = 0.025 70, fo- 
cussing on one of the two K points (fco = 2n/(3a)). (Dashed 
lines: two innermost bands for W — 40 x \fHa.) Here we 
use hopping matrix elements 71 = 73 = 0.12 70, but set 
74 = t' = 0. The inset shows the transverse charge dis- 
tribution of two edge states with energies indicated by the 
markers. Lower panel: The same but including hopping terms 
74 = 0.05 70, t' = 0.0770. 



the two layers, which is obtained from V g by including 
the screening in the layers 4 We assume that V is uni- 
form across the system, even near the boundaries, which 
in reality will be enforced to a good approximation by 
the proximity of the metallic gate close to the system. 
At the boundary, we impose standard hard-wall bound- 
ary conditions, which amount to setting the wavefunction 
to zero on lattice sites that lie outside the ribbon. 

It is also conventional to include the direct coupling 
73 between the non-dimerised atoms in both layers^ Ig- 
noring for the moment other hopping terms, a typical 
dispersion relation of the edge states in a gated ribbon 
is shown in Fig. [T] (middle panel). Of the four bands of 
edge states, two are massive, with a gap A v — 2\V\ (the 
same as for the bulk states), while the other two dis- 
play a linear dispersion, corresponding to two counter- 
propagating states which are well localized at opposite 
edges (see inset). The same scenario is replicated at the 
other K point in the graphene Brillouin zone, but with 
the edges interchanged. Therefore, each edge supports 
two counter propagating state, one from each K point. 
Notably, to this level of approximation, the propagation 
velocities on both edges are the same. 

The middle and lower panels of Fig. [T] shows how the 
dispersion changes when the two most prominent addi- 
tional hopping elements are taken into account: the in- 
terlayer hopping 74 ss 0.05 70 between adjacent dimerized 
and non-dimerizcd carbon atoms, and the intralayer hop- 
ping t' ss O.O770 between next-nearest neighbors. Both 
terms break the particle-hole symmetry, which then dis- 
criminates the different edges: the terminating atoms live 
on different layers and therefore posses a different on-site 
potential. In the dispersion relation, this takes the effect 
of an additional background velocity, which increases the 
velocity at one edge (the 'fast edge') while reducing the 
velocity at the other edge (the 'slow edge'). The depen- 
dence of these velocities on the gate potential is shown in 
the top panel of Fig. [5J The additional hoppings also af- 
fect the transverse localization of the edge states. This is 
shown in the bottom panel in terms of the inverse partici- 
pation (IPR) Iipr = Yl n l^™l 4 ' which is larger the better 
localized a state is (here, the sum is over the transverse 
direction of the nanoribbon). 

At this point it is important to mention that there 
are two types of zigzag edge termination in graphene 
bilayer^ the so called a and /3 terminations. When the 
on site potential V is zero some qualitative differences in 
the bandstructure can be found in both cases. However 
when V g leads to a nonvanishing V in both cases two of 
the edge bands become gapped (with a gap of the order 
of V) and the two other bands cross each other around 
the Fermi points as discussed before for the two types 
of edge terminations. Also in both cases the asymmetry 
between the velocities at different edges can be found, 
being the differences most prominent at the high energy 
bands. Thus we will not consider any distinction between 
the two terminations as long as we are only concerned in 
the physics at small energies around the Fermi points. 
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FIG. 2: Influence of 74 and t' on the edge-state velocities (top 
panel, vd = y/3yoci/2h) and localization on each isolated edge 
(Ah = 0) (bottom panel, expressed in terms of the inverse 
participation ratio -Zipr). 



At small energies, the hybridization of the edge states 
leads to an avoided crossing, where the small residual gap 
Ah vanishes exponentially with increasing width of the 
ribbon. A sizable hybridization gap exists only for 
small values of W (narrow ribbons) . For typical applied 
gate voltages and widths, each edge of the nanoribbon 
therefore behaves like a metallic one-dimensional system 
with two counter-propagating states, where each state is 
associated to one of the K points. In this scenario of 
low dimensionality, one should expect that interactions 
play an important role. Generally, however, interactions 
are felt most strongly when particles propagate slowly 
and are well confined. Thus, the concrete manifestations 
should depend on the distinct properties of the fast and 
slow edge states described above. This is what we will 
explore in the remainder of this paper. 



III. MANY BODY EFFECTS 

Within the tight-binding description, interactions can 
be incorporated via a Hubbard term U ni^n^i, where 
ni,a = c\ a Ci. a is the occupation operator on lattice site 
i with electrons of spin a =1, f- Assuming Ah = 0, each 
edge can be treated separately via a hamiltonian of non- 
chiral bosons: 



dx 



2K„ 



dx 



(1) 



where \i labels the charge (p) and spin (a) sectors, while 
Vfj, and are the velocity and the Luttinger parame- 
ters, respectively. Interaction events that do not lead to 
backscattering renormalizc these parameters in different 
ways, which yields the usual separation between the spin 
and charge degrees of freedom^: 
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where the upper sign applies to p, while the lower sign 
applies to a . 

Most of the interesting physics arises from backward 
scattering, which leads to a nonlinear Sine-Gordon term 
in the spin sector, 



-Hback — 
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27t 2 a 2 



dx cos(Vs<f> a ). 



(4) 



At half filling, which we need to consider to assess the 
size of any interaction-induced gap, this is accompanied 
by umklapp processes, which contribute an analogous 
term to the charge sector: 



umklapp 
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2n 2 a 2 



dx cos 



(5) 



All these terms are controlled by the forward and back- 
ward coupling constants which are obtained by writ- 
ing the Hubbard interaction term Ufii^rii^ in the basis 

9{x,y) ~ ^(y)e- lk ° x L(x) + p + (y)e ik °* R(x) (6) 

of the diagonalized one-particle hamiltonian around each 
K point (fcrj = 27r/ (3a)): 

9f = Ua 2 J dy\v±(y)\ 2 W ± {y)\ 2 , (7a) 

g b = Ua 2 [ dy\ V ±{y)\ 2 \^(y)\ 2 - (7b) 



The left-right symmetry of the system results in the ad- 
ditional constraint gt = = g = UalipR, which in- 
volves the IPR of the transverse wavefunction. As seen 
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above, as soon as 74 and t' are taken into account the 
IPR strongly depends on the layer-symmetry breaking 
potential V. Thus, the effective strength of interactions 
can be controlled by the applied gate voltage. 

Since the terms Q and (0 make the hamiltonian not 
exactly solvable, we assess their consequences using stan- 
dard renormalization group (RG) arguments j^2i In first 
order in gb, the RG equations for y M 



and K,, are 



dK p 
dl 



2 • 



dl 



= (2-2K fl )y ll , 



(8) 



where I is a logarithmic renormalization scale. The scal- 
ing behavior of y p depends on the value of K^. For 
Kfi > 1, y p decreases when the scaling parameter / in- 
creases, which renders the interaction term irrelevant. 
For K p < 1, however, the interactions are relevant; 



the 

intermediate case K p = 1 represents a quantum critical 



point. 

Because of the constraint <?/ = 
parameter ^ can be written as 



9b = g, the Luttinger 
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For the case g > (repulsive interactions), to which we 
restrict our attention, K„ > 1 and the RG equations (|5]) 
imply that the spin sector will flow towards the nonin- 
teracting theory, with renormalized K* = 1 due to the 
emerging spin-rotation invariance. On the other hand, 
K p is smaller than one, and therefore the system may 
acquire a gap in the charge sector. For a more quantita- 
tive analysis, we introduce the new variable x p = 2K p — 2, 
so that Eq. takes the form 



dx p 
~dl 



= ~Vp 



dy P 

dl 



(10) 



In these new variables, the RG equations possess a first 



integral A 



b p,0 



yi o, whose value remains 



constant along the trajectories of the flow, and therefore 
can be calculated using the bare values x Pt o and y Pi o- 
Most importantly, none of the flow lines cross the line 
x p = y p , where A = 0. This line constitutes a separatrix 
between two regimes. For y p < x p , y p scales to zero and 
x p approaches a definite renormalized value x*. Phys- 
ically, this sector of the system is then well described 
by the noninteracting hamiltonian (fTJ) with renormalized 
K*. For y p > x p , however, y p — > 00 and x p — ¥ — 00 flow 
towards strong coupling, where the RG equations are no 
longer valid. The charge sector then acquires a gap A p . 

For bare value \x P fl\/y Pi o -C 1, this gap can 
be estimated using the selfconsistent harmonic 
approximation, 2 - which exploits that the gap scales 



with ! as A, 



v p A. 



This leads to 



, = A e , with A 

1 

( 4 (aA)£ ) 2 ' where we equate the real- 
space a cutoff with the lattice constant. The quantity A 
is an ultraviolet cutoff in the momentum space. Since 
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FIG. 3: Charge gap A p as a function of V , for selected values 
of U, assuming 74 = t' = (upper panel), as well as 74 = 
O.O570, t — O.O770 (middle panel: slow edge; lower panel: 
fast edge). Note the different scales of the vertical axis. 



the edge dispersion relation is linear only for energies 
less than Ay ss 2V, A ~ J^-, which delivers 



A p = 2V 2 ~ 2KP 



(11) 



However, this estimate is only valid when \x p o\ <C y p ,o *C 
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1. For \x p .o\ < y Pi o <C 1, the gap must be obtained 
by terminating the flow at y p ~ 1, where the first-order 
expansion breaks down, giving 

A p = 2T/ v /TT^e- 1 ^ = 2V»)e~^- (12) 

In order to determine which of these regimes applies 
to the bilayer graphene edge states, we have computed 
the ratio x/y for some values of V and U. Consistently, 
\Xpfi\/lJpfi ^5 1: which means that the gap should be ob- 
tained from Eq. ()12|) . where v and g = UalipR follow 
from the results of Fig. [2] 

The size of the resulting gap for different scenarios is 
shown in Fig. [3] In all cases, the dependence of A p on V 
is non-monotonic. Ignoring next-nearest-neighbor hop- 
ping (left panel), the gap at both edges is identical, but 
takes a sizeable value only for unrealistically large values 
of the Hubbard parameter. When 74 and t' are taken 
into account, the gap at the fast edge (right panel) is 
further suppressed. However, the gap at the slow edge 
(middle panel) is dramatically increased, in particular in 
the region where the propagation velocity becomes small. 

IV. CONCLUSIONS 

We have studied the effect of next-nearest-neighbor 
hopping and interactions on the edge states in gated 
bilayer graphene nanoribbons with zigzag termination. 
The additional hopping results in the formation of a slow 
and a fast edge, where electrons propagate with different 
velocities. The small velocity at the slow edge assists the 
formation of an interaction- induced gap A p . Whether 
a sizeable gap can be achieved at realistic interaction 



strengths would be best decided by an experiment which 
exploits that umklapp processes are suppressed when the 
temperature is raised. The physics behind any experi- 
mentally observed gap can therefore be probed via the 
temperature dependence of the conductivity (measured 
at the edge via side contacts) , which for T ^> A p should 
scale as o oc T^~ ak p^- The experimental set up can be 
a four-terminal geometry which probes each edge sepa- 
rately (Ah = and thus the two edges are decoupled 
for a sufficiently wide nanoribbon). Because each edge 
state possesses different a different Fermi velocity K p 
is thus different as it is discussed in the text allowing 
to measure independently each conductivity. This alge- 
braic dependence is distinct from the exponential behav- 
ior of the conductivity when the gap originates in the 
single-particle spectrum (such as the hybridization gap 
Ah due to the finite width of the ribbon). The simul- 
taneous presence of the hybridization gap and interac- 
tions could be modelled by a bosonized interaction term 
# gap = -^ r sin(v / 8'/) p )cos(v / 8'/) CT ) r&2i which couples 
the spin- and charge sectors. Moreover, in this situation, 
right- and left-moving states around the same K point 
become coupled because of their finite overlap. The fact 
that these states propagate with different velocities adds 
interesting complications, which, however, go beyond the 
scope of the present work. 



V. ADKNOWLEDGEMENTS 

We gratefully acknowledge discussions with F. Guinea, 
M. A. H. Vozmediano and G. Leon Suros and support by 
the EC via Grant No. MEXT-CT-2005-023778. 



1 T. Ohta et al, Science 313, 951 (2006). 

2 E. McCann and V. I. Fal'ko, Phys. Rev. Lett. 96, 086805 
(2006). 

3 E. McCann, Phys. Rev. B 74, 161403(R) (2006). 

4 E. V. Castro et al, Phys. Rev. Lett. 99, 216802 (2007). 

5 H. Min, G. Borghi, M. Polini, and A. H. MacDonald, Phys. 
Rev. B 77, 041407(R) (2008). 

6 L. Jiao et al., Nature 458, 877 (2009). 

7 L. C. Campos et al., Nano Letters 9, 2600 (2009). 

8 M. Y. Han et al., Phys. Rev. Lett. 98, 206805 (2007). 

9 E. V. Castro et al., Phys. Rev. Lett. 100, 026802 (2008). 

10 J. Rhim and K. Moon, J. Phys.: Condens. Matter 20, 
365202 (2008). 

11 B. Sahu, H. Min, A. H. MacDonald, and S. Banerjee, Phys. 
Rev. B 78, 045404 (2008). 

12 M. P. Lima, A. Fazzio, and A. J. R. da Silva, Phys. Rev. 
B 79, 153401 (2009). 

13 A. Rycerz, J. Tworzydlo, and C. W. J. Beenakker, Nat. 
Phys. 3, 172 (2007). 



14 L. M. Malard et al, Phys. Rev. B 76, 201401 (2007). 

15 Z. Q. Li et al, Phys. R ev. Lett. 102, 037403 (2009). 

16 A. B. Kuzmenko et al.. larXiv:0908.0672l 

17 K. Sasaki, S. Murakami, and R. Saito, App. Phys. Lett. 
88, (2006). 

18 C. Kane, L. Balents, and P. A. Fisher, Phys. Rev. Lett. 
79, 5086 (1997). 

19 L. Brey and H. A. Fertig, Phys. Rev. B 75, 125434 (2007). 

20 T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford Science Publications, Oxford, 2003). 

21 A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, 
Bosonization and Strongly Correlated Systems (Cambridge 
University Press, Cambridge, 1998). 

22 T. Giamarchi, Phys. Rev. B 44, 2905 (1991). 

23 C. M. Varma and A. Zawadowski, Phys. Rev. B 32, 7399 
(1985). 

24 L. Balents and M. P. A. Fisher, Phys. Rev. B 53, 12133 
(1996). 



